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We study quantum mechanically the self-similar black hole formation by 
collapsing scalar field and find the wave functions that give the correct semi- 
classical limit. In contrast to classical theory, the wave functions for the black 
hole formation even in the supercritical case have not only incoming flux but 
also outgoing flux. From this result we compute the rate for the black hole 
formation. In the subcritical case our result agrees with the semiclassical tun- 
neling rate. Furthermore, we show how to recover the classical evolution of 
black hole formation from the wave function by defining the Hamilton-Jacobi 
characteristic function as W = Mining. We find that the quantum corrected 
apparent horizon deviates from the classical value only slightly without any 
qualitative change even in the critical case. 
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I. INTRODUCTION 



Matters can strongly interact gravitationally to form black holes. The Schwarzschild, 
Reissner-Nordtrom, Kerr and Kerr-Newman black holes are supposed to describe the end 
state of such gravitationally collapsing objects and are characterized by mass, charge, and 
angular momentum. Since these black holes are stationary, it would be physically interesting 
to understand the dynamical process for black hole formation that leads to the stationary 
black hole as the end state. Such an attempt toward understanding the black hole formation 
process was first done by Christodoulou who analytically studied a massless scalar field to 
prove that it either forms a black hole for the strong self-gravitational interaction or disperses 
for the weak self-gravitational interaction The gravitational collapse of the massless 
scalar field was confirmed in an early numerical simulation 0]. The most important results 
of gravitational collapse was first discovered by Choptuik who found through numerical 
investigation that an initial Gaussian packet for the massless scalar field collapses self- 
similarly from spatial infinity either to form a black hole or to disperse back into spatial 
infinity depending on whether the parameter characterizing the wave packet is above or 
below a critical value, and that the black hole mass exhibits a power-law scaling behavior 
whose critical exponent is independent of the initial data ||. Since then critical phenomena 
of gravitational collapse in the spherically symmetric geometry have been found numerically 
in the perfect fluid complex scalar SU(2) Einstein- Yang-Mills []10| |, non-linear a 

axion-dilaton [I2~,r3||, and gravitational wave models . 



The simplest model for the gravitational collapse is a massless scalar field evolving self- 
similarly in the spherical-symmetric geometry, and has been studied analytically Q and 
numerically @|§. The analytic solutions [15] for that model together with the self-similarity 



has been used to study the critical behavior [16-18]. Most of these works treat the classical 



aspects of gravitational collapse. Since the massless scalar field model predicts the type-II 
critical behavior, the mass of black hole near the critical parameter can be arbitrarily small. 
So one may ask whether quantum effects play a significant role in the system. Hence it would 
be physically interesting and important to study the quantum gravitational collapse. As a 
quantum mechanical treatment of gravitational collapse, the wave function of a supercritical 
parameter for a quantum black hole was expressed in terms those of the subcritical param- 
eters and thus explains the black hole decay quantum mechanically jn5J. Recently, we have 
also studied such quantum effects that, though classically the collapsing scalar field below 
the critical parameter does not form a black hole, quantum mechanically a black hole can 
be formed through tunneling process [p0|] . However, the interpretation of the wave functions 



in the supercritical and subcritical cases has not completely been settled yet fll9| , |20|] , and a 
question has been raised quite recently how quantum effects can change the physics such as 
the black hole mass or critical exponent near the critical value |H| . 



It is the purpose of this paper to study quantum mechanically the black hole formation 
in the massless scalar field model and to investigate how quantum effects modify the clas- 
sical picture of gravitational collapse. We use the ADM formulation to quantize the model 
and find analytically the black hole wave function in terms of the confluent hypergeometric 
function. The quantum model resembles in many respects a quantum FRW cosmological 
model minimally coupled to a massless scalar field. The standard interpretation of quantum 
mechanics and quantum cosmology enables us to calculate explicitly the rate for black hole 
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formation. We put forth a quantum criterion on the black hole formation, whose semiclas- 
sical approximation agrees with the classical one. Furthermore, by using the characteristic 
function defined by W = him In ip we are able to compute the quantum effects to the classical 
apparent horizon and to the critical exponent. 

The organization of this paper is as follows. In Sec. II, using the ADM formulation we 
derive the Wheeler-DeWitt equation for black hole formation by a self-similarly collapsing 
scalar field. The Wheeler-DeWitt equation has an SU(2) group structure and the wave 
functions are found in terms of the confluent hypergeometric function. In Sec. Ill, we 
find the wave function for the black hole formation, which consists of both the incoming 
and outgoing components. We calculate the incident and transmitted fluxes, in terms of 
which the rate for black hole formation is found. It is found that the rate for black hole 
formation gives the correct semiclassical limits in both the supercritical and subcritical 
cases. In Sec. IV we use the new method for recovering the classical solution from our wave 
function introduced by us recently in Ref. The key idea is to modify the Wheeler- 

DeWitt equation to have a nonvanishing energy E rather than the usual vanishing value. 
Then we define the Hamilton- Jacobi characteristic function as W = hlmlnip, from which 
we get the Hamilton- Jacobi evolution of classical solution. Finally we let the energy E 
vanish. Following this procedure we regain the semiclassical limit of our quantum solutions 
and find the position of the apparent horizon. We have used the steepest descent method 
to determine the position of apparent horizon and found good agreement with numerical 
evaluation. Even in the critical case (co = 1) there appears no qualitative change from the 
classical results. In the last section, we discuss the salient points of our work. 



II. CANONICAL QUANTIZATION 



The model to be studied in this paper is the spherical-symmetric geometry minimally 
coupled to a massless scalar field. Then the Hilbert-Einstein action 



S 



1 



lQn Jm 



d xy/—g 



i^-2(V0) 2 



+ 



8nG JdM 



d 3 xKVh, 



(1) 



where G is the Newton constant, reduces to the (1 + l)-action 



S sp h = - I d a;V-7 r 



^{^(7) + ^((Vr) 2 + l)}-2(V# 



(2) 



where 7 a & is the metric in the remaining two-dimensional manifold. The (1 + l)-metric for 
the scalar field collapse has the form 

{2) ds 2 = -2e 2a dudv (3) 

where u and v are null coordinates, and a is a function of u and v. To find the continuously 
self-similar solutions, r = —u^pzy = \/—uvy, <fi(z) and cr(z), which depend only on z = 



-2t 



we introduce a diagonal gauge 



u = —ue 



v = cue 



(4) 
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in terms of which the (1 + l)-metric takes the form 



^ds 2 



,2(7 



-2uo 2 dr 2 + 2duj' 



(5) 



The classical solutions with vanishing a are found in Refs. [|T^,|T^], and the corresponding 
spacetime has the geometry M. = R 1,1 x S 2 . 

In order to canonically quantize the system we adopt the ADM formulation and introduce 
a lapse function N(t) to write the (1 + l)-metric as 



-2N 2 (r)uj 2 dr 2 + 2du 2 . 
After taking into account the boundary terms, we obtain the actionF] 



(6) 



Ssph = J d(-^- —)dr 



1 



y 2 + 



2N 



hy 2 



2m 2 P N 



^ 2 -N{\y 2 



w d \2 2h ( d a 2-, 
l-^-y) +— ( u —d>) } 



mi 



(7) 



where overdots denote derivative with respect to r. Note that K = -jr\ — fr-£r has the 
dimension of action h and is proportional to the area factor measured in the Planck length 
lp, and plays the role of a cut-off parameter of the model. Hence ^ is a dimensionless 
parameter. The canonical momenta for y and are given by 



K 



HKy 2 
m 2 P N ( 



(8) 



Then Eq. (^) can be rewritten in the ADM formulation as 



S. 



where 



sph 



nj 1 2 , m P J2 

U = ~2K n y + 2hKf n * 



dx 



(9) 



(10) 



is the Hamiltonian. The N acts as a Lagrange multiplier, so one gets the Hamiltonian 
constraint 



n = o. 



(11) 



According to the Dirac quantization method, the constraint (|TT|) becomes a quantum con- 
straint on the wave function 



x In this paper we use a specific unit system c = 1 keeping explicitly the Planck constant h and 
the gravitational constant G, in which mp = \p^,lp = VhG, and mp = The variables u and 
r have the dimension of lp and the time-like variable r and y become dimensionless. 
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7^,7^0,7^)1*) = 0. (12) 
We further impose the constraints from the self-similarity 

gw-o-gm. da. 

Finally we obtain the Wheeler-DeWitt equation for the quantum black hole formation 



h 2 d 2 hm 2 p d 2 / 
2Kdy~ 2 ~~ 2Ky 2 8<P 2 ~ v 1 2 



y- 



#(</,</>) = 0, 



(14) 



which is the same equation that was used to calculate the tunneling rate for black hole 



formation in the subcritical case | 2"0fl . 

The wave function can be factorized into the scalar and gravitational field parts 

Kc n 



^{y,4>) = exp(±i . 4>)i>{y)- 
v n 1 mp J 



(15) 



Here the wave function for the scalar field was chosen to yield the classical momentum 
Kff, = h ^?N = ±Kc , where c is a dimensionless parameter. Then Eq. (|T4] ) reduces to the 
gravitational field equation 



2K dy 1 



K 1 - — 

2 



2r 



1>(y) = 0. 



(16) 



It is worthy to note that Eq. fllBD , as one-dimensional quantum system, describes an inverted 
Calogero model with energy K and has the group structure SU(2) with the basis |2c 



2 



n 2 K 2 c 2 t TTyy + yity 



(17) 



Due to the group structure we are able to find the solutions to Eq. ([H]) in terms of the 

(18) 



confluent hypergeometric function p4 
ipi(y) = Di 

where 



iK , 



(^-y 2 y-M(a.,b.^y 2 ) 



1 

^ = 4 




1 

a - = 2 




b_ = 1 


"5* 



with 



• > /r \' 2 



= (i^ - — 



(19) 



(20) 
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The other independent solution is given by 



ipn(y) = D 



•K 



K 



exp( 



2 h 



(21) 



However, the wave function ([Zl]) is a linear combination of ipi an d another solution 



in 



•K 



K 



cxp( 



2 h 



n 



(22) 



where 



1 i „ 



(23) 



In the above equations D's are constants, and /i's, a's and 6's are dimensionless parameters. 
Hence Eqs. (0), ( PI] ) and ( f2"2"D are dimensionless functions. 



III. WAVE FUNCTION FOR BLACK HOLE FORMATION 



In the classical context, the massless scalar field imploding self-similar ly from spatial 
infinity with momentum above the critical value collapses to form a black hole without 
leaving any remnant and with momentum below the critical value it reflects back to spatial 
infinity |TTJ. However, in the quantum context, we may follow the analogy to the scattering 
problem of quantum mechanics. We prescribe the boundary condition for the wave functions 
for the black hole formation such that they should be incident from spatial infinity and some 
part of them be reflected by the potential barrier back to spatial infinity but the remaining 
part be transmitted toward the black hole singularity inside their apparent horizons. 

We now wish to calculate the quantum mechanical rate for black hole formation. It is 
not difficult to see that Eq. ( 



has the asymptotic form p4| at spatial infinity 



BH 



My) = Dj 



(K A 



2 r y 



+ 



r(fe-) 
r( Q _; 



K , 



iK y 2 



e 2 « 



, (24) 



where Dj = Dj(—i)^-. To show that Eq. ( |T8| ) satisfies indeed the boundary condition for 
the black hole formation, we note that the first term describes the incoming component and 
the second term describes the outgoing (reflected) component at spatial infinity, (y ^> 1). 
Near y = 0, Eq. (|T8|) has the form |24|] 



if, 



BH 



h 



exp 



(25) 



so it has only the incoming flux toward the black hole singularity as will be shown later. 
We calculate the incoming (incident) flux at spatial infinity 



6 



Im 



.71",, 



- Qcoshf h (Q-K) 
1 11 (HK) 1 / 2 sinhfQ 



£(2Q+K) 



(26) 



where we used the relations |r(l + ia;)| 2 
hand, the outgoing (reflected) flux at spatial infinity is similarly found, 



-f^and|r(i + ,x) 



COsh WX 



241. On the other 



Jref. 



\Dif 



Qcosh^(Q + ijQ ^ 
(^)V2 S i n h |Q 



K 



(27) 



The rate for black hole formation is the ratio of the transmitted flux to the incident flux. 
From the flux conservation applied to Eqs. ( |2^D and fl2"7p we obtain the (transmission) rate 
for the black hole formation 



Jtran. - Jref. 

Jin. Jin. 

cosh ^(Q + if) _ f( 
cosh§(Q-^) e " 



1 - 



(28) 



Eq. (p8|) determines the probability for the black hole formation for all values of Cq. 

In order to understand the classical limit of this solution, let us consider the large K 
limit. 

From the mass of 



K always enters Eqs. (|18|), fl24|), and fl28|) in the combination o: 
limit is equivalent to the small h limit, which is the semiclassical limit 
black hole formed in the supercritical case (T7| 



m 



h V2 



K 



we can express y as 



M AH = -f- -r AH 



K 



m 2 P 



rrip (UJ C \ 



Cq ^ Dip ' 



(29) 



(30) 



where u; c denotes the cutoff of the uj variable. In the supercritical case the influx from spatial 
infinity falls into the black hole singularity and the black hole mass increases infinitely for 
an infinite duration. Hence the large K limit is the large black hole mass limit. In this limit 
the rate for black hole formation should be calculated separately for Co > 1 and Co < 1. For 
c > 1 Eq. (|28|) becomes asymptotically 



Jtr 



-f^(co-l) 



Ji 



(31) 



Therefore, the rate for black hole formation becomes unity for Co > 1, which implies that a 
black hole is formed from the collapsing scalar field, as in the classical case IE 
In the case of c < 1, we carefully rewrite Eq. (|^) as 



Jtr 



e -^(K+Q) ginh |Q 

cosh^(fr-Q) ' 



(32) 
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The asymptotic value of Eq. fl32|) 



Jtran - = e ~l K ( 1 -^)^ (33) 

Jin. 



is the tunneling rate for quantum black hole formation in the subcritical case ||20|| . The 
classical result is obtained in the very large limit of K, in which the scalar field with cq > 1 
evolves supercritically to form a black hole and with Cq < 1, it bounces back without forming 
a black hole. 

A few comments are in order. Firstly, the flux conservation is valid at spatial infinity 
and the singularity. To show the conservation we compute the transmission flux from spatial 
infinity toward the singularity 

Jtran. Jin. Jref. 

= Wl 2 ^. < 34 > 

which coincides with the flux near y = obtained by a direct computation using Eq. fl25|), as 
expected. Secondly, since Eq. (|14|) is linear and Co is just a separation parameter, one may 
construct a more general solution for the black hole formation by superposing the solutions 
([24]) with different c 

^bh(vA) = f dc F(c )ij B H(y)exp(±i Jj C ° 4>), (35) 
J v ft 1 mp J 

where F(co) is an arbitrary weighting factor. By recalling that (/i__ — a J) from Eq. (|l^) and 
(/!_ — a* + ) from Eq. (^) are independent of Cq, one still has the factorized asymptotic form 
forEq. (k§ 



Here, -F(±)(0), which are obtained by integrating over cq, may represent wave packets of 

± ._ Ken , 



plane waves, e R m p , in Eq. ([T5j) . 

Finally, we explain the reason why the other wave functions for Eq. ([14]) are irrelevant 



for the black hole formation. The wave function (|2^) has incident components from both 
spatial infinity and the black hole singularity to the potential barrier, so it does not satisfy 
the boundary condition for the black hole formation. Another wave function is obtained by 
taking a linear combination of Eqs. (|18|) and ( p2|) which describes the outgoing component 
alone at spatial infinity, 

-^--Sgl^SH(4,T--" ; -^. (37) 



r(a_)v n 

Since the net flux is the difference between those of the incoming and outgoing components, 
it is easy to see that the one-parameter (8) family of wave functions, 



S 



,-,W 



- 1 



]_ gi7r(a_ — a+) 



r(a-)r(a») 
r(a+)r(a;) 



7V 



(38) 



has the same rate as Eq. (|28|). However, the branches and ^f>jy of the wave function fl3"8"|) 
have the incident components both from spatial infinity and from the black hole singularity 
to the potential barrier, respectively. We also note that the wave function 4>} v , the complex 



conjugate of Eq. (|37D , is also a solution to the linear real equation (jT4|) and is purely 



incident from spatial infinity. However, near y — 0, consists of two branches of wave 
functions: ipj having an outgoing flux and ip* ni having an incoming flux. Thus %p* IV describes 
the incident waves from both spatial infinity and the black hole singularity to the potential 
barrier. Hence these wave functions, ip* IV and xjjg(9 ^ 0), do not satisfy in a strict sense the 
boundary condition for the black hole formation. 

IV. SEMICLASSICAL LIMIT AND APPARENT HORIZON 



In order to regain the semiclassical picture from our quantum wave function we use our 



recently proposed method [ 22] . For this we first consider the Wheeler-DeWitt equation with 
nonvanishing energy parameter E given by 



h 2 d 2 hm 2 P d 2 



$ E {y,<f>) = E^ E (y, 



(39) 



2Kdy 2 2Ky 2 dcj) 2 

where E has the dimension of K. The solution corresponding to the black hole formation is 



Mv) = D 



%K , 



i^y 2 T M{a E ^\y 2 ) 



(40) 



where /i_ and 6_ are the same as Eq. fll9D , while 

1 i 



CLE 



-2H( Q + K + E 



(41) 



Note that the parameter E appears indirectly only through cle- The Hamilton- Jacobi char- 
acteristic function Wsiy) is given as 



W E (y) = hlm\nip E (y) 



From Eq. 



1 1 K K 

= —-Ky 2 — —Q ln(—y 2 ) + hImM(ciE, b-, i^-y 2 ) + constant. 

Li It lb 

we recover the evolution equation of the gravitational collapse as 



(42) 



E=0 



h 



_d_ 
~dE 



K 



Im In M(a E , b-,i—y 2 ) 



E=0 



(43) 



where (3 is a constant to be determined by an initial condition. 
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Before turning to the evolution equation ( [43|) in the region where quantum effects are 
important, we compare it with the classical equation of motion in the regions where classical 
effects are dominant. The classical solution of the evolution equation 



is given by 



^ + 1-^-^ = 
n 2 2y 2 ' 



r + (3 = ± j dy 



y 2 + f-i 



±\\n(y 2 - l + v / y 4_ 22/ 2 + Cj 



Firstly, at spatial infinity (y — > oo) we have the asymptotic expression 



(44) 



(45) 



t + (3 = ±- In y 2 + constant. 



(46) 



We compare this classical evolution with the quantum one. We compute the characteristic 
function for the incoming and out-going components, separately, using Eq. 



t + (3 = him 



— ^ <9 „ r . , .K 9 , 
2MTa M(a > b ->V 



a=a— ,y—*oo 



(47) 



= ±^y 2 , 

where the upper (lower) sign comes from the outgoing (incoming) component in Eq. (p|). 
Hence in the asymptotic region the Hamilton- Jacobi solution and the classical limit give the 
same evolution of the collapsing process. 

Secondly, near the singularity (y — > 0) the allowed motion is toward the singularity y = 0. 
So the classical solution is given by the lower sign in Eq. ( f45"D and has the asymptotic form 



r + /3 = --m 



co-l + f 1 



1 

Co 



h constant. 

2c 



(48) 



On the other hand, the wave function ( ^0|) depends on E only through the confluent hy- 
pergeometric function M(a^, 6_, i^-), and the semiclassical limit of the quantum solution 
is 



r + f3 = h 



8_y 

dE 

2 



K 

Im In M(cie, b_,i—y 2 
n 



y 

2c 



(49) 



Thus, in the semiclassical limit (h — > 0, which is equivalent to K — > oo) the quantum wave 
function has the correct limit to the Hamilton- Jacobi solution both in the asymptotic region 
and near the singularity. 
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In the intermediate region between the asymptotic region and the singularity it is most 
interesting to see how the apparent horizon and black hole mass are affected by quantum 
effects. We note that the position of the apparent horizon is given by the trapped surface 



(Vr) = 0, which translates to y — y ■ 
determined from the equation of motion, 



0. In the classical case the apparent horizon is 



dW c 
dy 



A 



K 2 (y 2 + C \ 



Ky = -Ky, 



(50) 



which has the solution 



v 2 ah 



(51) 



On the other hand, in the quantum case we use the quantum counterpart Eq. ( f42"D to Eq. 
(0), 



dW F 



dy 



E=0 



Ky = —Ky. 



Since W\e= 



hlmlnip, Eq. fl52|) becomes 
d 



— Ky = h— (im In ijj 

dy y 



Q 



-Ky h him 

y 



dy 



(52) 



(53) 



which can be rewritten as 



Q_ 

2K 



hy 2 Re 



aM(a + l,b+l,z) 
b M(a,b,z) 



(54) 



In the semiclassical limit {K/h — > oo) we can evaluate the apparent horizon {y\n) using the 
steepest descent method (see the Appendix for details). We find the apparent horizon up to 
the first order as 



Vah 



c 2 h 2 3 

— H 

2 A 2 



64 + 80 h 3 

cq 2 _ cq 4 _|_ n( ' 
8 l K 3 ' 



(55) 



where the first term in the right hand side is the classical result and the second term is 
the first order quantum correction. We checked this with numerical calculation for large 
K/h(= 40). They are 



Vah 



= 1.9616(c = 2); 0.5000(c = 1); 0.2569(c = -=) 

numerical \/1 



vah 



Eq{55) = 1.99 94(c = 2); 0.5001(c = 1); 0.2652(c = -j=). 



(56) 



Even in the critical case (cq = 1) the apparent horizon deviates only slightly from the classical 
one, and the first order correction is quite good. As far as the apparent horizon is concerned 
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the critical case seems to show no special behavior compared to the super and subcritical 
cases. For small K/H there is a large quantum correction, and the approximation (|55j) is not 
valid. For example we calculated numerically when K/H = 1, and the deviation from the 
classical value is substantial as 



- y 2 

numerical 



classical = -°- 412 ( c o = 2 ); 0.127(c = 1); 0.128(c = -^=). (57) 



However, there is no qualitative change in the position of the horizon. 



V. DISCUSSION 

In this paper we have studied the black hole formation by a self-similarly collapsing 
massless scalar field. The analytic wave functions for the black hole formation were found in 
terms of the confluent hypergeometric function. By evaluating the incoming and outgoing 
flux at spatial infinity we were able to compute the rate for quantum black hole formation, 
which agrees with our previous result in the subcritical case. To find other quantum effects 
we used the characteristic function defined by the imaginary part of the wave function and 
recovered the evolution of black hole formation in time with quantum effects taken into 
account. We now compare our results with those from the classical solution, and the wave 
functions by Tomimatsu, and discuss the physical implications. 



Firstly, in the classical case |T7| the evolution of the gravitational collapse is governed 
by Eq. and the solution is given by Eq. (|45|) . In the supercritical case the incoming 
component, the lower sign of Eq. P5|), falls into the singularity and any part of it is not 
reflected back to spatial infinity. On the other hand, in the quantum case the wave function 
(|I%|) for black hole formation has both the incoming and outgoing components (^). In 
particular, in the supercritical case it has a small fraction of outgoing component. In the 
subcritical case the black hole formation is allowed by quantum mechanical tunneling. 

Secondly, among all the wave functions (|T8|), (f2~l|) and (|22| ) for the Wheeler-DeWitt 
equation, Eq. ([Tj|) has the desired semiclassical limit for black hole formation. On the other 
hand, Eq. (|3T| ) represents a wave function with outgoing component alone. Though the 
physical meaning of this wave function is not clear, it may be interpreted as an amplitude 
for the black hole decay or white hole creation, which requires a further study. 

Thirdly, as mentioned in Sec. IV, the apparent horizon in the quantum case deviates 
from the classical one only slightly. Thus there appears no qualitative change from the 
classical to the quantum treatment of the apparent horizon. 

Finally, we compare our wave functions with those by Tomimatsu who found the wave 
functions for the quantum black hole formation and decay [19| . The difference is that he used 
a different coordinate system from ours and quantized both the classical constraint equa- 
tion and the Hamiltonian from the reduced action, whereas we adopted directly the ADM 
formulation, obtained the Hamiltonian constraint and quantized the constraint according 
to the Dirac quantization method. In contrast to his wave functions diverging at the origin 
in both the supercritical and subcritical cases, our wave functions are always regular at the 
origin. The regularity of wave function for black hole seems to suggest that quantum gravity 
effects may cure some singularity problems in classical gravity. 
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APPENDIX: USEFUL FORMULAS FOR THE CONFLUENT 
HYPERGEOMETRIC FUNCTIONS 



The asymptotic form of the confluent hypergeometric function |2~5[ used in this paper is 

M(a, 6, z) e iwa z- 



r(6) 



r(6-a) 



1 + 



e z z a-b 

T(a) 



(Al) 



37T 



where — | < arg z < 



To obtain the characteristic function valid at intermediate region we use the integral 

\b-a-l 



formula for the confluent hypergeometric function [p^| 



r(&-a)r(a) 



r(6) 
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dte zt t 
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ty 



(A2) 



which is valid for Re b > Re a > 0. In this paper z = ij^y 2 and Re b = 1 > Re a = | > 0. 
In order to evaluate M(a, b, z) in the large K limit we evaluate it with the steepest descent 
method. For this we define 



iK 

—fit) =tz + (a-l) ln(t) + (b-a-1) ln(l - t). 
n 



We find the relevant root of % = as 



2 - 9 ?/ 2W2 



_ c + y 2 - (y 4 + c 2 - 2y 
2y 2 

where c = (cq 2 — j^) 1 ^ 2 — ij?- We evaluate the integral 



(A3) 



(A4) 



dt exp(i—f(t) 
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After straightforward calculations we obtain 
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*2 o _ 64 , 80 



CO 



n 2 



K 2 



4Cn 2 



(A5) 



(A6) 



which is in good agreement with the numerical calculations in the large K limit as shown 
at the end of Sec. IV. 
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